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Abstract 

The relationship between the lower limit on the nuclear stability lifetime as 
derived from the non disappearance of 'stable' nuclei (Tf^ ^ 5.4 x 10^^ yr), 
and the lower limit thus implied on the oscillation time {Tnn) of a possibly 
underlying neutron-antineutron oscillation process, is clarified by studying 
the time evolution of the nuclear decay within a simple model which respects 
unitarity. The order-of-magnitude result Tnn ~ 

2(rd/r^)i/2 > 2 X 10^ sec, 

where Tn is a typical n nuclear annihilation width, agrees as expected with 
the limit on Tnn established by several detailed nuclear physics calculations, 
but sharply disagreeing by 15 orders of magnitude with a claim published 
recently in Phys. Rev. CRAP. 
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I. INTRODUCTION 



The stability of nuclei, as determined by looking for proton decay places a lower 
limit on the 'disappearance' lifetime > (5.4 ± 1.1) x 10^^ yr, and sets a lower limit also 
on the lifetime of other hypothetical processes such as neutron-antineutron (nn) oscillations 
in free space. The lower limit on the free-space oscillation time Tnn which has emerged from 
several quantitative nuclear physics calculations using n and n optical potentials 0-^] is 
approximately given by 

r„^^2(Trf/r^)^/2>(2.1±0.2)xl0^sec , (1) 

where ~ 100 MeV is a typical n-nuclear annihilation width. This slightly exceeds the 
limit 

Tnn > 0.9 X 10^ sec (2) 

set in the ILL-Grenoble reactor experiment |0. A new experiment planned at ORNL will 
hopefully improve this limit by two orders of magnitude ||^. 
Eq. may be rewritten as 

showing that, since TnTnn ^ 1 owing to the huge lifetime distinction between strong interac- 
tions (lifetime ~ T^^) and superweak interactions (lifetime ~ Tnn)-, the nuclear disappearance 
lifetime induced by nn oscillations is many orders of magnitude longer than the nn oscilla- 
tion lifetime in free space. This is equivalent to the common statement that nn oscillations 
in matter undergo severe suppression, since the neutron and the antineutron feel nuclear 
potentials which are extremely different from each other, and the mass degeneracy which 
allows such pure oscillations between a free n and a free n is thus removed in the nuclear 
medium. 

Nazaruk [mTU| has raised objections to the use of nuclear physics potential models, claim- 



ing that nn oscillations are not suppressed at all in the nuclear medium. Thus, in the first 
paper he obtained ~ Tnn, resulting in the limit Tnn > 10^^ yr which advances the lower 
limit (|l|) by about 30 orders of magnitude! Dover, Gal, Richard |11] and Krivoruchenko |12 



subsequently pinpointed errors in that paper and reestablished the generally accepted lower 
limit (|1|). In the second, more recent paper |T^, Nazaruk 'rederived' the ~ Tnn result. 



but argued that since the nuclear decay is non exponential, a more careful consideration of 
the nuclear stability limit translates into Tnn > 10^^ yr, which advances the lower limit (|I]) 
by 'only' 15 orders of magnitude. Such far-reaching claims, when published by a first-rate 
journal, should not go unanswered. 



In this note I wish to expand the arguments outlined by Dover et al. ||IT|, in order to 
show in detail how Eq. (|I]) for Tnn is obtained. In Sec. |I| we consider the nn mass matrix in 
n-annihilating matter, in order to derive in the most economical way the eigen-lifetimes in 



a toy seesaw model. In Sec. |T| we study in detail the time evolution of the nn 'oscillating' 
system in matter. In the concluding Sec. |IV| we have tried to identify the error in Nazaruk's 
logic. 
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II. OSCILLATIONS THE SEESAW MECHANISM 



A common approach in problems of oscillations, be it for CP symmetry or in neutrino 
mass generation problems, is to diagonalize the mass matrix in order to find out the physical 
states. For nn oscillations induced by a coupling e, the simplest form of the in-medium mass 
matrix is 



m e 
e m — 



(4) 



where m is the joint value of the n and n masses, and ~ 100 MeV is the annihilation 
width of the n in the nuclear medium. For F^ = 0, i.e. in free space, the two eigenmasses 
are m±e, differing only by a tiny e from each other. The reality of these eigenmasses gives 
rise to a purely oscillatory behavior n n, with lifetime Tnn = In the nuclear medium, 
when Tfi ^ 0, the eigenmasses of (H) are given by 



16e 



m 



eigen ^ 4^^™ i 4^-^?! ( 1 -p2 ) ' (^) 



F? 

n 



which to leading order in £^/F? assume the values 



mn = m- -^^Y? ■ ■ ■ ' ^n = m- -iVn + . . . , (6) 

where the dots stand for higher order terms. The degeneracy of the real parts of meigen can 
easily be removed by introducing n-nucleus and ra-nucleus real optical potentials which differ 
substantially from each other. In discussing decay modes, however, only Im meigen matters; 
adding real potentials does not change the splitting between, and the order of magnitude 
of the two values for Im meigen in Eq- (H)- Therefore, for the sake of simplicity, we ignore 
these real potentials. Identifying the eigenwidth 7 with —2 Im meigen, one obtains from Eq. 
(P) two eigenwidths: 

Ae 

In = , 'yn = Tn . (7) 

J- n 

Intuitively, since = F^ is the nuclear annihilation width of an antineutron, 7„ must stand 
for the decay rate of the neutron (in units where ^ = 1). It is clear that this decay rate un- 
dergoes a huge suppression factor, 4e/F^ -C 1, with respect to the free-space oscillation rate 
e. The huge disparity between 7^ and is a good demonstration of the seesaw mechanism 
encountered in the discussion of neutrino mass generation problems; note that the product 
Inln = Ae^ is of the same order of magnitude as in free space, which is a necessary condition 
for the seesaw mechanism. These qualitative arguments need to be explored quantitatively 
by studying the time evolution of nn 'oscillations' in the nuclear medium. This is done in 
the next section. 



HI. EXPLICIT TIME EVOLUTION 



For simplicity, following Ref. |11 



we write the time dependent coupled Schrodinger 
equations for zero momentum neutron and antineutron in nuclear matter {Ti = 1): 
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(8) 



No nuclear (real) potentials Un and Un appear in the present discussion, which is focussed 
on the interplay and competition between the free-space oscillation rate e = and the n 
decay rate Tn due to the strong-interaction n nuclear annihilation, and their effect will be 
briefly discussed later on. The coupled first order equations (|^) give rise to the following 
second order differential equation for each one of ipnyi^n- 



(9) 



Seeking eigensolutions of the form tpj = exp{iujjt/2), the eigenfrequencies uj satisfy a 
quadratic equation 







(10) 



yielding two solutions: 



iTf 



^ 2 



'/2 



^r,i -4- + ... 



(11) 




As' L 

1 + 4— + ... 



(12) 



where the expansion in even power of the extremely small quantity /T\ ^ 1 is indicated. 
These eigenfrequencies are related to the eigenwidths in Eq. (^ by = 27. The eigenfre- 
quency ojn gives rise to exponential decay of l^/'p with a rate F^,, to leading order. Dover et 
al. [|ll|] noted that ojn should be discarded since it corresponds to the rate of disappearance 
of an antineutron when its oscillation coupling e to the neutron is neglected. It was further 
noted that the eigenfrequency u;„ gives rise to a rate of disappearance 'i{e/Vn)e which was 
then interpreted as that for a neutron, due to its oscillation coupling e to the antineutron. 
This latter rate is equivalent to the estimate given in Eq. (^ for T^. We note that since 
— 45^, the huge disparity between the frequencies Un and Un is a good demonstration 



of the seesaw effect encountered in the discussion of scales of masses in high energy physics. 

Here we wish to proceed in greater detail and rigor, using the eigenfrequencies (|TT], |12|), 
to construct the wavefunctions ipj with the appropriate temporal boundary conditions for the 
neutron (j = n) and for the antineutron (j = n). In considering the nuclear disappearance 
lifetime T^, these boundary conditions are 



tlJn(t = 0) = 1 , tlJnit = 0) = . 
The corresponding linear combinations of the eigensolutions are then given by 



^n(t) 



UOf. 



exp {iuJnt/2) 



- — exp {iuJnt/2) 



(13) 



(14) 
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2e 

i^n{t) = (exp {iuJnt/2) - exp {iuJnt/2)) . (15) 

Note that since ujn/ujn ~ 4£:^/r? -C 1, ipnif) is dominated at all times by the eigensolution 
with eigenfrequency a;„. It is also clear, by inspecting Eq. ([TT|), that the n wavefunction 
which evolves from zero at t = always remains very small with respect to the neutron 
wavefunction: 

|^.(t)/^„(t)r~^, (16) 

n 

this approximate relationship becoming exact for times t ^ F^^ when the second exponent 
on the right-hand side of Eqs. ([T^ , |15[) may be safely dropped out. For such 'long' times, 
the time dependence of both ip^ and ipn is given by 

I^Wl'~exp(^-4^t^ , (17) 

so that the rate of disappearance F^^ may be read off this exponential decay: 

4^2 

r. = T^ . (18) 

J- n 

Recalling that e = t~^, we obtain for the disappearance lifetime 

Td = F^ = — (FfjT^fj) Tnn ; (l"^) 

which agrees precisely with Eq. (||) of Sec. |I[ 

The minimal model of Eq. (P) for the time evolution of ipn and ipn may be extended to 
include real potentials f/„ and Un, and spatial structure: 

idti/jn = -^i'n + Un{r)'tpn + £^n , (20) 

2m 



idti'n = -T^^n + {Un{r) " iW {r)) + , (21) 
2m 

with the same t = boundary conditions specified in Eq. (|T^. In order to obtain the decay 
rate of the neutron, we multiply Eq. (pOD for ip^ by ■?/'*, and the corresponding equation for 
-0* by ipn, subtract from each other and integrate over space. The result is 

-dtj \iPn? dh = 2ejlm (rn^n) d\ . (22) 

This decay rate vanishes for £ — > 0, but the precise power of e by which it occurs depends 
on the ratio {i/Jn/i^nl- Acting similarly on Eq. ( ^T]) for ijjn, the decay rate of the antineutron 
is obtained: 

-dtj \^nfd'r = J 2W\iJnfd\-2e J Im (V^* rfV . (23) 
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We note that the 2e J terms in Eqs. 0) appear with opposite signs; what's lost out of 
the neutron intensity due to the oscillation coupling e is precisely gained by the n intensity, 
and vice versa, on top of the n annihilation decay rate. The total decay rate is obtained by 
adding up Eqs. (H, H): 



-dtj {\^nf + m')d^r = j2W\^nfd'r, (24) 

showing that if the antineutron did not annihilate in the nuclear medium {W = 0), there 
would be no loss of intensity from the combined n, n space. A rough estimate of |'?/'fiP/|V'nP 



may be made by inspecting Eq. (^) and noticing that, for the boundary conditions Eq. 
(P!B|), the evolved {Un — iW)ipn term on the right-hand side of Eq. ( PT| ) should be of the same 
order of magnitude as the source term eipn which generates it. Since Un and W ~ Tn/2 are 
all of the same order of magnitude, we obtain 



(25) 



which agrees with Eq. (^6|). The properly normalized decay rate 7 is then given, using Eq. 
(0), by 



-3, /(!</'„ 








/( 




2 






(Pr 



7 = ,/| , , , ,2\ = . /, , ,2 . ,2\ ' (26) 

which is approximately of order 







r 








d^r 



7 ~ 2Vr 



f- , (27) 

-L r. 



n 



agreeing with the order of magnitude of F^, Eq. ([181) . We note that Eq. (|26D expresses the 
loss of n and n intensities to the unspecified final nuclear debris products. 

IV. DISCUSSION 

We outlined in the previous sections, using a simple and transparent potential model, how 
the free-space nn oscillation period gets tremendously prolonged in the nuclear medium, 
by a factor TnTnn/4:, to yield the corresponding nuclear decay lifetime of Eq. ([T9|) . More 
detailed calculations which treat the nuclear medium as a dynamical entity, confirm 
this order-of-magnitude estimate of T^. It is then perplexing to encounter occasionally claims 
that nn oscillations are not suppressed in the nuclear medium. Below we counter the most 
recent claim of this sort made by Nazaruk |llO[| . 



Nazaruk argued that potential models involve double counting by allowing for n-nucleus 
elastic and inelastic scattering, and that this is manifested within such models by a calculated 
nuclear decay probability which is linear in the time t instead of the (correctly anticipated) 
t"^ dependence for 'short' oscillation times. Our response is that the only nuclear property 
of the antineutron included in the present potential model is its nuclear decay width Tn, 
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unrelated in this model to any underlying n scattering process. Furthermore, we do recover 
a leading dependence by expanding Eqs. (14, 15) in powers of time t: 



l-\^IJn\' = eH^-]:e^Tnt' + ... , (28) 



\iJnf = eh'~^6^Tnt^ + ... , (29) 

so that the neutron depletion rate would appear to be e, as in free space — unaffected 
by the nuclear medium, in particular by the n annihilation width Tn, and agreeing with 
Nazaruk's thesis. However, this expansion is useful only for extremely short times, such 
that Tfit <C 1, before the n annihilation can have its act and exercise its toll. For relevant 
times which satisfy Tnt 3> 1, the expansion (^, ^9]) becomes useless. As argued here in Sec. 



pn] , however, for such (still short) times, the Un exponent in Eqs. (|T4|, |1^) may safely be 
neglected, resulting in the following exponential decays 

{t)f ^ exp (-4^^) , (t)l' ^ ^exp (-4^^) , (30) 

with a rate of disappearance = Ae'^/Tn (Eq. (|18])) which is directly read off this expo- 
nential decay. 

The trouble with Nazaruk's arguments may be traced as follows. Using S-matrix manip- 
ulations, for t ^ F=^^ he reaches the expression 

w{t) ^ eh^Wfiit) (31) 

for the 'probability' w{t) of nuclear decay in terms of the n decay probability Wn{t) (see Eqs. 
( P!T| , |I^) of Ref. |T^). He then assumes that, since over this time all the antineutrons got 
already annihilated, Wn{t) ~ 1. This is wrong: the n decay probability is a conditional one, 
depending on how many antineutrons were there in the first place with respect to neutrons. 
Unitarity prevents the accumulation of too many n, and as seen from Eq. ( pO[ ) 



(f,i « t « T,,/e') ^ ^ . (32) 



Hence 



w it) ^ ^t' , (33) 

n 

so that the neutron disappearance rate is of order 2e^/Ffj, in agreement with Eq. (^T]) and 
with our result Eq. (|T^ for F^. 

To summarize, we have confirmed the well known lower limit obtained on nn oscillations 
in free space from the stability of nuclei, as given by the estimate of Eq. (|l]) which states 
that Tnn ^ 2 X 10^ sec. This estimate agrees to within better than a factor of two with the 
series of quantitative calculations by Alberico et al. and Dover et al. 0. There is good 
reason then to plan reactor experiments |^ which promise to push up this lower limit by 
perhaps as much as two orders of magnitude. 



This research was partially supported by the Israel Science Foundation 



6 



REFERENCES 



[1] M. Takita et al. (Kamiokande collaboration), Phys. Rev. D 34, 902 (1986). 
[2] Ch. Berger et al. (Frejus collaboration), Phys. Lett. B 240, 237 (1989). 
[3] K.G. Chetyrkin, M.V. Kazarnovsky, V.A. Kuzmin, and M.E. Shaposhnikov, Phys. Lett. 
B 99, 358 (1981). 

[4] W.M. Alberico, A. Bottino, and A. Mohnari, Phys. Lett. B 114, 266 (1982); W.M. 

Alberico, J. Bernabeu, A. Bottino, and A. Molinari, Nucl. Phys. A429, 445 (1984); 

W.M. Alberico, A. De Pace, and M. Pignone, Nucl. Phys. A523, 488 (1991). 
[5] C.B. Dover, A. Gal, and J.M. Richard, Phys. Rev. D 27, 1090 (1988); Phys. Rev. C 31, 

1423 (1985); Nucl. Instr. Meth. A 284, 13 (1989). 
[6] J. Huefner and B.Z. Kopeliovich, Mod. Phys. Lett. A 13, 2385 (1998). 
[7] M. Baldo-Ceolin et a/., Z. Phys. C 63, 409 (1994). 

[8] Yu.A. Kamyshkov, in Intersections between Particle and Nuclear Physics, edited by 

T.W. Donnelly, AIP Conf. Proc. 412 (AIP, New York, 1997), p. 335. 
[9] V.I. Nazaruk, Phys. Lett. B 337, 328 (1994). 
[10] V.I. Nazaruk, Phys. Rev. C 58, R1884 (1998). 

[11] C.B. Dover, A. Gal, and J.M. Richard, Phys. Lett. B 344, 433 (1995). 
[12] M.I. Krivoruchenko, Phys. At. Nucl. 59, 1972 (1996). 



7 



